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Spiral waves are ubiquitous in two-dimensional systems of chemical or biological oscillators coupled 
locally by diffusion. At the center of such spirals is a phase singularity, a topological defect where 
the oscillator amplitude drops to zero. But if the coupling is nonlocal, a new kind of spiral can 
occur, with a circular core consisting of desynchronized oscillators running at full amplitude. Here 
we provide the first analytical description of such a spiral wave chimera, and use perturbation theory 
to calculate its rotation speed and the size of its incoherent core. 
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Large systems of coupled limit-cycle oscillators have 
been used to model diverse phenomena in physics, chem- 
istry, and biology [H, H|. Examples range from spiral 
waves in heart muscle 3[ and brain tissue to syn- 
chronization of nerve cells [5[ , pendulums [f| , fireflies [7| , 
and Josephson junctions [8[. 

Most of the early theoretical work on these systems 
assumed either local coupling (through nearest-neighbor 
or diffusive interactions) or global coupling (through 
infinite-range interactions, corresponding to a mean-field 
approximation). In the past few years, however, several 
researchers have begun to explore other types of connec- 
tivity. One line of research investigates what happens if 
the coupling is neither local nor global but somewhere in 
between, as occurs naturally in certain neural, chemical, 
and biochemical systems 0, [HI [H| • 

Around 2002, Kuramoto and his colleagues discovered 
a spatiotemporal pattern that appears to be unique to 
such nonlocally coupled systems [T3, E3 ■ They showed 
numerically that systems of identical oscillators with 
symmetrical coupling could self-organize into a state in 
which some oscillators were mutually synchronized while 
others remained desynchronized. An exact solution for 
this "chimera state" was later obtained for phase oscil- 
lators arranged on a ring (l3l . [l4j , and its dynamics and 
bifurcations were further clarified by studying two inter- 
acting populations of oscillators O [l(| [TtJ . 

The most surprising result emer ged from simulations 
of two-dimensional arrays [H, Il9l. l2Cl|. Kuramoto and 
coworkers [TH, [2(| [2l| showed that two-dimensional 
chimeras are spiral waves with a phase-randomized core 
of desynchronized oscillators surrounded by phase-locked 
oscillators in the spiral arms. Nothing like this had ever 
been seen before; in ordinary reaction-diffusion systems, 
the smoothing effects of diffusion would prevent the spa- 
tial discontinuities implied by an incoherent core. 

In this Letter we present the first analytical solution for 



a spiral wave chimera. Our analysis predicts the existence 
of the incoherent core and yields a simple formula for its 
radius as well as for the rotation rate of the spiral arms. 
The model we consider is 



dt 



= lj- I G{\x-x'\)sin[^(x,t)-^(x',t) 



y]dx 
(1) 

where <f>(x,t) is the phase of the oscillator at position 
x £ R 2 at time t, uj is its natural frequency, and a is a 
phase lag. The kernel G is a normalized Gaussian: 

3 -|x-x'| 2 



G(|x-x'| 



(2) 



Note that the width of G determines a length scale for the 
system. Varying a can be thought of as choosing different 
first harmonics of a general coupling function [22| . In 
contrast, u> is dynamically irrelevant, and can be set to 
zero without loss of generality. 

Figure QJ a) shows a spiral wave chimera for Eqs. {TJ, 
(|2"|) , computed numerically. Notice the phase-randomized 
core at its center, as highlighted in Fig.[]7b). The spatial 
structure of this pattern is governed by an integral equa- 
tion, originally derived by Shima and Kuramoto [l8| via 
the following self-consistency argument. Define a local 
mean field (a complex order parameter) by 



R(x,t)e 



iO(x,t) _ 
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so that (UJ) can be rewritten as 
d<f>(x,t) 



dt 



lu - R(x, t) sin [0(x, t) - Q(x,i) + a]. (4) 



Now seek statistically steady solutions of Eqs. ([I]), J2}. 
For this special class of solutions we can find a rotating 
coordinate frame (in state space, not physical space) in 
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FIG. 1: Snapshot of a spiral wave chimera on a square do- 
main of side length 15, using an array of 151 x 151 oscillators. 
Parameter values: a = 0.157T, uj — 0. Colors schematically 
indicate oscillation phases in panels (a)-(c), whereas in panel 
(d) they depict the amplitude i?(x) of the local order parame- 
ter according to the numerical values shown next to the color 
bar. (a): Spatial variation of sin<ji>. (b): Close-up of the in- 
coherent core in panel (a), (c): sin(0(x)) (phase of order 
parameter) (d): -R(x) (amplitude of order parameter). The 
amplitude drops to zero at the center of the core, consistent 
with the randomized phases of the oscillators there. 



which the phase-locked oscillators in the spiral arms ap- 
pear motionless, and the mean field becomes stationary 
across the entire pattern [l8|, . To define this rotating 
frame, let ip(x, t) = 0(x, t) — ilt and = 6 — ilt where 
f2 is as yet unknown. Then using ^ and ((3]) we find 



dt 



w-fi-#(x)sin[<0(x,t)-e(x)+a!] (5) 



where the time-independent mean field is given by 
i?(x)e i6(x) = / Gdx-x'De^'^dx'. 



(6) 



Oscillators for which \oj — Sl\ < i?(x) will reach a steady 
state, ip*(x), f° r which dip/dt = in Eq. ([5]), whereas 
those for which \u> — 0| > i?(x) will drift but dis- 
tribute themselves with a stationary density in the ip 
variable inversely proportional to their angular velocity 
dip/dt [3, [II]. Equation ([6|) can then be solved self- 
consistently by replacing e 1 ^ by e 1 ^ (for locked oscilla- 
tors) and by the mean of e 1 ^ (for drifting oscillators), 
where the mean is calculated using the density just men- 
tioned. The result is the nonlinear integral equation 

iJ( x ) e *eW = ie~ la j G(|x-x'|)e ie ( x ') 
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where A = w — fi. Figures UX C ) andQJd) show sin8(x) 
and i?(x) for a typical spiral wave chimera. 

Our new results concern the analytical solution of 
Eq. l[7]). The first step is to simplify it using the spi- 
ral wave ansatz of Cohen et al. [23|. Considering the 
appearance of the functions O and R in Figs. [TJc) and 
(d), it is natural to use the ansatz i?(x) = A(r) and 
O(x) = 6 + *S?(r), where (r, 9) are polar coordinates. Sub- 
stituting this into (0 one obtains 



A(r)e 



oo 

i< ■- / K (r, s)e i * (s) 
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(9) 



and I\ is the modified Bessel function of the first kind of 
order 1 . This ansatz has replaced the problem of finding 
the functions of two variables, i?(x) and O(x), by that of 
finding A(r) and ^(r), both functions of a single variable. 

Nevertheless, solving Eq. ([8]) is still a formidable task. 
To make progress, we use perturbation theory to ana- 
lyze (|8|) in the limit of small a. First observe that at 
a = 0, Eq. ((T|) is a gradient system in a coordinate 
frame rotating at speed lo; hence, in this frame all the 
attractors for fTJ) must be fixed points. Back in the orig- 
inal frame, these fixed points correspond to phase-locked 
states rotating at frequency f2 = u>. Therefore A = 
when a = [lj|, which motivates the series expansions 

A = Aia + 0(a 2 ) (10a) 
A(r) = Ao(r) + Ai(r)a + 0(a 2 ) (10b) 
*(r) = *oW + *i(r)a + 0(a 2 ) (10c) 



as a -> 0. Substituting (fTOa) - (fTUc]) 
lowest order, 



into (TST) we obtain, to 



i4o(r)[l + i* (r)]= / K(r, s)[l +W (s)} ds. (11) 
Jo 

Hence 

/•OO 

Mr) = / K(r,s)ds 
Jo 

= (V^/2)re- r3 /2[j o(r 2 /2) + Il(r 2 /2) ] (12) 

and 

nOO 

Ao(r)V (r)= K{r, s)^ Q {s) ds. (13) 
Jo 

Equation (fl~2")) shows that at leading order, the ampli- 
tude A (r) of the mean field increases with r and satis- 
fies Aq(0) = and linv-joo Aq{t) = 1, in agreement with 
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the behavior seen when a / (Fig. [ljd)). Furthermore, 
under the physically reasonable assumption that v I / o( r ) 
can be written as a power series in r, one can show that 
^oM = C is the only solution of (fTB")) : then, from the ro- 
tational invariance of the problem we can set x f'o( r ) = 0. 
This result indicates that the spiral arms approach radial 
straight lines as a — > 0, a fact we have also confirmed nu- 
merically. 

At 0(a) we obtain A\(r) = and 



4>M*i(r)= / K(r,s) 
Jo 



4>(«) 



ds- A (r) 



which, after defining /(r) = Ao(r - )$i(r), yields an inho- 
mogeneous Fredholm equation of the second kind: 



f(r) 



K(r,s) 



f(s)ds = Ax 



K(r, s) 



ds-A (r). 



Ao(s) 

(14) 

We have been unable to solve (jT4j) analytically, but 
extensive numerical investigations suggest that to sat- 
isfy /(0) = (which we know to be true, given that 
Aq(Q) — 0), Ai must be very close to 1. To show that 
this is plausible we note that for r* ^> 1, 



A (r*)«l and K(r*,s) 



,-(r*-s 



(15) 



so that if f(r) is slowly varying near r* , the left side 
of (fl4|) is approximately zero while the right side of (fl"4"l) 
is approximately Ai — 1. 

Recalling that it is those oscillators with A(r) < |A| 
which drift, we can see that to order a the radius of the 
incoherent core is p, where Ao(p) — |Ai|a. In the limit 
as a — > 0, A (p) — > A' o (0)p = (y/n/2)p, and hence, using 
Ai = 1, we have 



a + 0(a 2 ). 



(16) 



Thus, to order a we have two testable predictions regard- 
ing the behavior of the spiral wave chimera as a — ► 0: 

1. Spiral arms rotate at angular velocity f2 = u — a. 

2. Incoherent core radius is given by p = (2/^/n)a. 

Figure [2 shows measurements of both A = oj — and the 
radius of the incoherent core from a simulation of Eq. (fTJ), 
as a is varied. The data points agree very well with the 
above predictions, despite the use of a finite domain. 

Along with testing our analytical results against sim- 
ulations, we have also tested them against numerical so- 
lutions of the integral equation {8}. The results for a 
domain < r < 15 are shown in Fig. [3] as a is varied. 
Figure [J3a) compares Aq (r) , the only non-trivial zeroth 
order term in our series expansion, to the computed A(r) 
obtained for two different values of a. Boundary effects 
near r = 15 are clearly visible. For a — 0.017T, A(r) 
is indistinguishable from Ao(r) for r < 13. Figure [3Jb) 
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FIG. 2: (Color online) Comparison of results from simulation 
and perturbation theory, (a): Values of A from simulation 
(circles) compared with predicted values A — a (solid line), 
(b): Comparison of measured values of p (circles) with p = 
(2/y/Tr)a (solid line). The simulation was run on a square 
domain of side length 10, with 201 x 201 oscillators. 
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FIG. 3: (Color online) Numerical solutions of Eq. (|8]). (a): 
Air) for two different values of a, and Ao(r). (b): ty(r)/a 
for three different values of a. (c): Calculated values of A 
(circles), compared with the predicted result of A = a (solid 
line), (d): Radius of incoherent core p: solutions of A(p) = A 
(circles), compared with the predicted value p = (2/y / 7r)a 
(solid line). 



plots \&(r)/a for three different values of a; this choice of 
dependent variable is motivated by Eq. (|10c[) , which pre- 
dicts that ^(r) should be proportional to a, for a small. 
Figures [3fc) and rSfd) show A and p as functions of a, 
now found by solving ([5]) instead of simulating (jTJ) , in the 
same format as Fig. [2] The agreement with the data in 
Fig. O is excellent. 

The results above were derived using the Gaussian ker- 
nel In contrast, Shima and Kuramoto [IH used the 
modified Bessel function of the second kind, Kq, as a 
kernel. This coupling function arises naturally in a class 
of reaction-diffusion systems where active chemical el- 
ements are indirectly coupled by passive, fast diffusive 
agents Using a Gaussian allowed us to explicitly 

perform the integral in ^ and thus make analytical 
progress in the subsequent perturbation analysis. For 
Kq and other kernels one can evaluate the integral in © 
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FIG. 4: (Color online) Effects of different kernels. Panels (a) 
and (b) show A and p respectively for kernels 2nG(r) given 
by: (A) A'o(r); (B) 2e" r2 ; (C) e" r ; (D) 2H(r)H(l-r), where 
H is the Heaviside step function. 



numerically. In Fig. Q] we show A and p as a function 
of a for four different kernels, found by numerically solv- 
ing ©. Interestingly we see that while A does not seem 
to depend on the details of the kernel, p does. 

Several puzzles remain about the spiral wave 
21] . For example, we do not understand how it 



chimera 

bifurcates. When a is increased sufficiently far from zero, 
the chimera becomes numerically unobservable. Perhaps, 
for example, it loses stability though a subcritical Hopf 
bifurcation, or is destroyed through a saddle-node bifur- 
cation. Even when the chimera does exist, its core wan- 
ders erratically. This can be rationalized by observing 
that there is no preferred origin for the spiral in an in- 
finite domain, so it is neutrally stable with respect to 
translations. The finite number of incoherent oscillators 
act a source of effectively random "noise" which then jig- 
gles the spiral in these neutral directions. 

A second puzzle is that the spiral wave chimera 



apparently exists only if a is sufficiently close to 0, 
whereas the chimeras found in lower-dimensional sys- 
tems [H, [T3, [ID, [H, [TtJ exist only if a is sufficiently close 
to 7r/2. In physical terms, two-dimensional chimeras live 
near the gradient system limit, whereas one- and zero- 
dimensional chimeras live in the opposite regime, where 
the dynamics have an almost reversible or conservative 
character. At the moment we have no explanation for 
this difference. 

Although we have assumed identical oscillators in our 
analysis, this is not necessary. Spiral wave chimeras 
should persist in the presence of slight disorder. For ex- 
ample, Laing [ll| studied a variant of Eq. |T]) for non- 
identical oscillators and with a kernel G which was zero 
for |x — x'| larger than a particular radius. He used the 
recent ansatz of Ott and Antonsen [24], 25] to derive a dif- 
ferential equation for a variable closely related to the or- 
der parameter. For oscillators whose intrinsic frequencies 
uj are randomly drawn from a Lorentzian distribution, 
the steady states of this differential equation satisfy ([7]), 
where A = u — il + iD and D is the half-width-at-half- 
maximum of the Lorentzian distribution. The results in 
Fig. [3 agree qualitatively with those of Laing [l9( . 

The possibility of observing spiral wave chimeras in 
physical systems naturally arises. Nonlocal coupling of 
the form used here occurs in some networks of neu- 
rons [1,0, so such systems may be the best candidates. 
We leave the experimental observation of chimera states 
as a challenge to others. 

Research supported in part by NSF grant CCF- 
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